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I. INTRODUCTION 

One of the most critical challenges for spintronic de¬ 
vices, as the often mentioned spin-field-effect-transistor 
due to Datta and Das^, lies in the control of the carrier 
spin lifetime. The latter is limited by the spin relaxation 
and dephasing processes in semiconductors. The predom¬ 
inant mechanism of the spin relaxation in such devices 
is of Dyakonov-Perel type^. To extend the application 
of spintronic devices to the non-ballistic/diffusive regime 
with spin-independent scattering, it is of particular inter¬ 
est to find conditions for the electrons/holes in the semi¬ 
conductors which result in symmetries that correspond, 
according to Noether’s theorem^, to the conservation of 
spin. These symmetries enable to detect long-lived or 
even persistent spin states, i.e., states which do not relax 
in time. In structurally confined two-dimensional elec¬ 
tron gases (2DEG) such persistent solutions have been 
predicted for a special interplay between the linear-in- 
momentum Dresselhaus"^ and Bychkov-Rashba®’® spin- 
orbit coupling (SOC) by J. Schliemann et al, Ref. 7, and 
were extended by B.A. Bernevig et al., Ref. 8. These 
special states have been later confirmed by means of op¬ 
tical experiments®’The first type of SOC appears in 
semiconductors with broken inversion symmetry in the 
crystal structure (bulk inversion asymmetry (BIA)), the 
second one occurs when a structure inversion asymmetry 
(SIA) in consequence of an asymmetric confining poten¬ 
tial in the semiconductor heterostructure is present. 

In contrast to electron systems, the SOC is distinct and 
much more complex for holes although the underlying 
fundamental mechanism, described by the Dirac equa¬ 
tion, is the same. Since the conduction band, for most 
semiconductors, is an s-type energy band and the valence 
band is of p-type, the qualitative variation comes from 
the different total angular momentum, which is j = 3/2 
in the valence band, giving rise to heavy and light holes 
(HH, LH) and the split-off holes. As the mixing of HH 
and LH strongly influences the SOC, the reduction of 
dimensionality, like the 2D hole gases (2DHC) in semi¬ 
conductor heterostructures, has an immediate effect. 


This is due to the fact that the size quantization causes 
an energy separation between HH and LH states even at 
a vanishing in-plane Bloch wave vector which affects the 
strength of the HH-LH mixing and thus the SOC at finite 
k vectors. As a consequence, the magnitude of the SOC, 
especially the prefactor for Rashba SOC, depends sen¬ 
sitively on the confinement, as will be discussed in this 
paper. Conversely, the Rashba SOC in the lowest con¬ 
duction band is hardly affected by the size quantization. 
This is due to the s-type character of the energy band and 
also since the SOC is mainly determined by the energy 
gaps between the bulk bands. These energy gaps, how¬ 
ever, do not differ significantly if adding a confinement. 
For a proper description of the hole system, compared to 
electron systems, a substantially higher number of SOC 
terms is needed and requires approximations for an ana¬ 
lytical investigation at an early stage. 

Also, internal or external strain can yield significant 
consequences to the hole band structure. The reason is 
that it introduces additional couplings between the HH 
and LH bands, whereas - assuming a semiconductor with 
a direct bandgap - the conduction band is only indirectly 
affected due to the interaction with the strain-altered va¬ 
lence band.^®“^^ Moreover, the cubic crystal structure of 
the semiconductor has an imprint on the symmetry of 
the hole spectrum as can be seen in the warping of the 
Fermi contours and these always follow the strained crys¬ 
tal symmetry. 

Nonetheless, despite their complexity hole systems of¬ 
fer opportunities not available in electron systems and are 
particularly interesting for practical device applications 
for several reasons. First, the large effective mass m* of 
holes compared to conduction band electrons diminishes 
the kinetic term such that contributions from SOC be¬ 
come more important. Second, the p-wave character of 
the HH and LH states reduces the hyperfine interaction 
of the carrier spin with the nuclei. This allows in prin¬ 
ciple for long spin relaxation/dephasing times.An¬ 
other important aspect is the strength of the SOC in hole 
systems which can reach several meV in the splitting as, 
e.g., shown in CaAs/AlCaAs heterostructures^®’^®. All 
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these features of p-type systems facilitate a very effective 
manipulation of carrier spins and, hence, motivate fur¬ 
ther studies of hole gases in semiconductors as the one 
presented here. Moreover, since in a 2DEG the conserved 
spin quantities are always limited by a fc-cubic Dressel- 
haus contribution the question arises whether this is also 
the case for the 2D hole systems. 

An appealing continuation of the findings on spin¬ 
preserving symmetries in electron systems is the analysis 
of persistent spin states in hole systems as done recently 
in Refs. 20 and 21. However, these publications presup¬ 
pose materials with strongly restricted and unusual band 
structures. Following Ref. 21, a strainless sample with 
both Rashba and Dresselhaus SOC allows for the exis¬ 
tence of a persistent spin helix (PSH) in a 2DHG only 
in the case of a vanishing Luttinger parameter 73 with 

7 1 > 0 and 72 > 0.^^ Most of the semiconductors can 
only be properly described using a band model where 

72 < 73 ^^^^'‘(App. A), though. In Ref. 20 the PSH was 
found in the presence of finite strain and Rashba SOG 
where the condition for the Luttinger parameters is re¬ 
stricted by a different, however, also unusual condition 
72 = - 73 ^®- 

Another approximation which is often applied is to 
drop all invariants in the bulk Hamiltonian which lead 
to BIA and have relatively small expansion coefficients. 
This procedure is justified in bulk systems. In 2DHGs, 
this approach leads to a model Hamiltonian with both 
Rashba and Dresselhaus SOC being essentially cubic 
in momentum^^’^®, in contrast to 2DEGs where the 
dominant Rashba term is linear and the Dresselhaus 
term linear and cubic in momentum.However, recent 
publications^^’^® which are related to the seminal paper 
by Rashba and Sherman, Ref. 11, show that the relevance 
of the linear Dresselhaus terms in 2DHG has been un¬ 
derestimated: the abovementioned estimations are thus 
questionable and fail at least for the standard compound 
GaAs. 

In this paper, we present conditions for a wider range 
of semiconductors including strain, linear and cubic Dres¬ 
selhaus and Rashba SOC under which conserved spin 
quantities can be found. 

This paper is organized as follows. In the next section, we 
derive an effective two-dimensional heavy/light hole like 
Hamiltonian including strain, Dresselhaus and Rashba 
SOC. In Sec. HI, we derive the conditions for the exis¬ 
tence of a conserved spin quantity. Thereby we discuss 
its realizability and apply our findings to a prominent 
compound, InSb. Finally, we summarize our results. 


II. THE MODEL 

The aim of our investigation is to identify the appro¬ 
priate interplay between BIA (Dresselhaus SOC), a con¬ 
fining potential H(r) (build-in and/or external) causing 
SIA, and strain (either externally imposed using, e.g., 
the piezoelectric effect or induced by the epitaxial growth 


process) which gives rise to a conserved spin quantity in 
the hole system. To find analytic conditions we derive 
an effective HH/LH like 2x2 model, depending on the 
character of the up-most valence band. Starting point is 
a 4 X 4 model which is derived form the extended Kane 
model and includes the Luttinger Hamiltonian^®. 


A. Effective 4x4 Hole Hamiltonian 


Hereafter, we choose the coordinates to be x || [100], 
y |j [010] and x || [001]. We use the Luttinger parameters 
7 i, the bare electron mass mo and elementary charge e > 
0, the electric field £z^ the total angular momentum J for 
j = 3/2 and the symmetrized anticommutator {A, H} = 
[AB + BA) /2. 

The applied model which we use as a starting point for 
the investigation is an effective 4x4 Hamiltonian given 

by 

H = Hl + Hbia + Hs + V. (1) 


The first term represents the Luttinger Hamiltonian for 
HI-V semiconductors 
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-I- c.p. 


^73 


} + c.p.] 


( 2 ) 


where c.p. denotes the cyclic permutation of the pre¬ 
ceding indices. The second term 'Hbia accounts for the 
Dresselhaus SOC and is decomposed by the theory of 
invariants as^^ 
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"Wbia = ^ Ck {kx{Jx, Jy — J^} + c.p.) 

+ blf {{kx,kl - + c.p.) 

+ blf'"{{ka;,kl - kl}Jl A c.p.) 

+ blf'^iik^^kl + kl}{J^, Jl - J/} + c.p.) 

+ j2} + c.p.). (3) 

It includes a fc-linear term proportional Ck which 
is usually rather smalk®’^^ and thus often not 
considered^^’^^’^®’^®. An exemplary comparison of the 
remaining cubic contributions with the extended Kane 
model in App. B 2 shows that in a bulk system the term 
proportional to 64 /®” is the most relevant one. How¬ 
ever, in 2D size quantization leads to additional linear 
terms that may become significant in certain parameter 
regimes^^’^^’®®. We note that there are discrepancies in 
the perturbative determination of the coefficients 
as outlined in App. B 2. 

Furthermore, the effect of strain is described by the 
Bir-Pikus strain Hamiltonian 


"Hs = 5] 


^ ^ Ji d ^ ^ Cij {Jii J j} | • (4) 

i, 3^i 
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We assume the strain, with being the symmetric strain 
tensor, to be uniaxial in-plane or biaxial in-plane. As a 
consequence, an in-plane strain with the tensor compo¬ 
nents Cxx, ^yyi ^xy yields Only one extra strain compo¬ 
nent, Czz- Thus, we set txz = Cyz = 0 in Eq. (4). One 
should stress that slightly different notations occur in 
the literature for the deformation potentials a, b and d. 
Some relations between different definitions are listed in 
App. D 1. Here, we follow Ref. 20. Thus, the potentials b 
and d are positive whereas a can be positive or negative. 

Eventually, the potential V = 4^4-14 includes the con¬ 
fining potential Vc{z) and the external potential Ve{z). 
The latter causes structure inversion asymmetry (SIA) 
which induces Rashba SOC. We assume the potential V 
to depend only on the z-coordinate, which is pointing 
in the growth direction [001] of the semiconductor het¬ 
erostructure. More explicitly. 


Ve{z) =l4x4 • e£zZ, 

and an infinite square well of width L with 


Vc{z) — 14x4 


0 for z G [0, L], 
oo otherwise. 


( 5 ) 

( 6 ) 


kx = ky = 0, and 'H' is treated as a perturbation with re¬ 
spect to the appropriate inverse splitting 1/Ahi, (1/Ahh) 
between a HH like and a LH(HH) like subband. The en¬ 
ergy splitting Ahi is due to both the spatial confinement 
in the [001] direction and the imposed strain. 

According to the confinement, the eigenstates of Ho 
are given by \j,mj) \n), the product of the eigenstates 
of the total angular momentum J with j = 3/2, = 

±1/2 for LH and rrij = ±3/2 for HH and the subband 
index of z-quantization n € N*. The eigenfunctions of 
the quantum well in position space are given by (zjn) = 
\J2!L s\a{z n-K!L) which lead to the matrix elements of 
the ki and z operators given by 


{n\kz\l) 

{n\kl\l) 

{n\z\l) 


2inl ((-!)'+” - 

L in? - P) 



1 ) 


(1 — dnl), 


4niL((-l)'+"-l) 
. 7i-2(Z2-n2)^ 


for n = I, 
otherwise. 


(9) 

( 10 ) 

( 11 ) 


The eigenenergies of Ho (twofold degenerate) according 
to the sub-spaces are given by 


Notice that the discontinuity of the potential V/ can re¬ 
sult in non-hermitian matrix elements of k^ if not taken 
care. This problem can be resolved by a regularization 
procedure as shown in Ref. ? . Furthermore, contribu¬ 
tions due to boundary effects which result from the pres¬ 
ence of heterointerfaces are assumed to be small. Such an 
interface can allow for additional HH-LH mixing. 

A possible alternation of spin relaxation due to interface 
effects will be discussed elsewhere. 


B. Effective 2x2 Model for the First Subband 

In the following, we choose the basis states in such a 
way that the upper left block represents the HH the lower 
right block the LH subspace: 

nj _ ( ’^HH 'Hhh-lh \ 

\'Hlh-hh Him ) ' 

The confinement in z-direction allows for a further simpli¬ 
fication of the model to an effective 2x2 Hamiltonian us¬ 
ing quasi-degenerate perturbation theory (Lowdin’s par¬ 
titioning). The full Hamiltonian H is separated into two 
parts 

H=Ho+H', ( 8 ) 

according to App. B 1. The partition is, in general, not 
uniquely defined and different ways of splitting H are 
possible. For the given system, a meaningful decomposi¬ 
tion is the one which allows a projection on the subspace 
of a particular HH or LH like subband. Here, we select 
Ho to contain the diagonal elements of'H l±’Hs±V/(z) at 


EhkH = - ^(71 - 272) ( kl ) r ?, ( 12 ) 

Ei^n{n) = - ^(71 ± 272 ) {kl) V? ± 4 (13) 

2777-0 

where 5 = b {exx + ^yy — corresponds to the di¬ 
agonal shear strain tensor and {kl) = {l\kl\l). For 
simplification we subtracted an overall constant energy 
shift of AF = — — (3/2)6ez2. It be¬ 

comes apparent that strain effects as well as the con¬ 
fining potential lift the degeneracy of HH and LH bands 

at fc|| = ^kx + kl = 0. 

In the following, the in-plane strain contribution will 
be rewritten, according to Ref. 20, defining an in-plane 
strain amplitude jd and orientation 9 as^^ 

:= J^{^ {^^^ ~ ^vv) + 2*dea:y)- (14) 

The effective Hamiltonian for the lowest subband is now 
derived by quasi-degenerate perturbation theory (see 
App. Bl). Thereby, the lowest subband can be either 
a HH like or a LH like subband, depending on the mag¬ 
nitude and type of strain. Up to third order in the energy 
splitting it yields 

3 

(15) 

P—0 

= + (16) 
P—0 

with the Pauli matrices (7^ and the H^^l where the 
superscript (p) indicates the order in the perturbation 
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(App. Bl). Additionally, we will neglect terms of the 
order oiO{k% 

We would like to mention some properties of the per¬ 
turbation theory first. Concerning the decomposition of 
the Hamiltonian, it should be stressed that if we instead 
choose Hq = Hh + Hs, that is, without the confinement 
potential, we can make use of the non-commutativity of 
the momentum operator kz and the position operator z to 
derive a finite Rashba coefficient as done in Refs. 20 and 
24. However in this case, the energy splitting Ahi cannot 
be a result of the subband quantization, but only have 
strain as an origin. This approximation may be justified 
if the strain splitting is much larger than the subband 
splitting. Yet, since we consider a quasi 2D system, the 
subband splitting is an essential effect and we are thus 
to choose the partitioning as described above. 

An important observation of the perturbation pre¬ 
sented in this paper is that a finite Rashba spin-orbit 
(SO) field, resulting from the coupling between different 
subbands, can only be obtained by third or higher order 
perturbation theory. This is due to the fact that the di¬ 
agonal elements of H', Eq. ( 8 ), and thus are not yet 
involved in the second order, according to Eq. (B7). In 
addition, it will be shown in the following that given a 
HH(LH) like ground state it is necessary to include, in 
addition to the first LH(HH) like subband, also the sec¬ 
ond HH(LH) like subband in the perturbation procedure 
to yield a finite contribution due to Rashba SOC. 

Concerning the significance of the various Dresselhaus 
contributions in Eq. (3) we mentioned earlier. Sec. H A, 
that in the bulk system keeping only the cubic term pro¬ 
portional is a good approximation. This is due 

to its large value compared to the other cubic terms. 
However, the size quantization causes an additional linear 


Dresselhaus contribution for HHs that would not appear 
if only the term proportional was considered. For 

the light holes, the situation is different as the term pro¬ 
portional 64 ) 1 ’®” already yields a linear term which clearly 
dominates over the remaining linear contributions. As a 
result, we will take into account the effect of the fc-linear 
terms generated in the first order perturbation by the 
terms proportional Ck, 642 ®’^ and 65 )'®’^ only in case of a 
HH like ground state and neglect them in case of a LH 
like ground state. The coefficient 652 ®’' yields only a small 
cubic term which can be disregarded. Based on its large 
value, for higher order BIA corrections we incorporate 
solely terms proportional to 644 ®’'. 

Depending on the nature of the strain, the splitting 
can lead to either a lowest HH like or LH like subband. If 
we, e.g., specify to the case of uniaxial compressive stress 
in [ 110 ] direction, we obtain d < 0 since exx = ^yy < 0 
and Ezz > 0.^® Therefore, the splitting between HH and 
LH is enhanced and the topmost subband is HH for 
arbitrary material since 7 ^ > 0. Consequently, we have 
to distinguish two cases: 

1. System with a Heavy-Hole Like Ground State 

Assuming the ground state being HH like and applying 
third order Lowdin perturbation theory, we get for the 
part of the Hamiltonian which is proportional to unity in 
spin space, according to Eq. (16), 

Kff,HH = 14ff,+ (Ahl,h2), (18) 

where we defined 


Ekin,± = — 


(71 ± 72 )fcj| ± ^ 

2mo " Ahi 


11,11 14 ^ 


8v8v 

'41 




— (^kx — ky) cos(2d) -I- 2kxky sin(26*) 
.73 


(19) 


and 


14ff.±(A) = 




3. 

-^bezz- 


{kl) 

2 mo 


(71 T 272) 


e£zL 


± 


3^4 


pH 


4m^Ahi ,11 


256L2e2£:| 

sTa 


( 20 ) 


The energy gaps are given by Ain,hm = EiM{n) — E'hh(w) and Ahn,im analogously. The effective vector field El due 
to Rashba and Dresselhaus SOC, which is modified by the presence of strain, yields 


rinn = 


(21) 
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with components given by 

^X,± = Ad.± kxkl (72 T 273) - klj2 + 0^13 




rj±kx + Ar,± [± (72 ± 273) klky =F 72^^ + /3^73 (0ky cos(26») + k^ sin(26»))] 


kly,± = Ad,± < ± klky (72 T 273) T kl-f2 + /3^73 


^ j sm(26») =F ky 1 ) cos(26») 


± r]±ky ± Ar,± [± (72 ± 273) kxkl T l 2 kl + /3^73 {ky sin(26») =F k^; cos(26»))] 


0^,± = 0. 


(22a) 


(22b) 

(22c) 


Here, the Dresselhaus coefficients r]±, Ad,± and the 
Rashba coefficient Ar^± are given by 


-!+ = - - 5 (1,’r 

,- = - btf" ikl) 


l8v8v 

^51 


Ad,± = ± 


Ar,+ — - 


2TOoAhi,ii 
128/i^e£z73 


'41 


'{kl}, 


Ar,- = 


97r2TO§Ahi.iiAhi.h2’ 

128S'^ef^73 

Q-TT^TO^All hlAll 12 ' 


) (kl) 


(23) 

(24) 

(25) 

(26) 
(27) 


The index (±) distinguishes the case of a system with a 
HH like ground state (+) from the one with a LH like 
ground state (—). 

Notice that given a vanishing 73 the dominant contri¬ 
bution due to Rashba SOC vanishes, too. Only a con¬ 
tribution as a result of the coupling to conduction bands 
is left, which is of higher than third order: In the bulk 
system, for most semiconductors, the dominant invariant 
in the extended Kane model^®, which is present due to 
SIA, is given by^"^ 


= rlf0{ky£. - k,£y)J, + C.p.). (28) 


If the bulk system is reduced to an effective two- 
dimensional system, the according counterpart of this 
term in an effective 2x2 Hamiltonian can be calculated 
using Lowdin perturbation theory as done above, keep¬ 
ing the factor r®)'®’' unchanged. However, as mentioned in 
Ref. 21, for a HH like ground state this resulting term is 
of higher order than the one given proportional to Ar (al¬ 
though represented by the same invariants). This can be 
understood by recalling the root of the coefficient r®)'®’^: 
It is the coupling between valence and conduction bands. 
In contrast, if a confinement is present, the contribution 
resulting from Rashba SOC in the effective HH system 
is dominated by the splitting between HH and LH like 
subbands. In the case of a LH like ground state an addi¬ 
tional /c-linear Rashba term proportional to r®)'®’^ appears 
already in third order perturbation theory. The angular 
momentum matrix is zero in the HH subspace but 


has finite matrix elements in the LH subspace. Since the 
prefactor r®)"®^ contains terms which are inversely pro¬ 
portional to the band gap^^, this contribution can be 
neglected since we assume the conduction bandgap to 
be much larger than the subband splitting. The contri¬ 
bution stemming from BIA has a different nature: The 
parameter 6®)'®’', which is connected with the invariant 
{{kx,kl — kl}Jx + c.p.) in the Kane model, is mainly 
defined through the valence band Tgi, and conduction 
band Tgc gap Eq. Thus, it is hardly affected by the sub¬ 
band quantization. Moreover, in contrast to the Rashba 
contribution, the corresponding Dresselhaus term in the 
confined system appears already in second order of the 
applied perturbation. Hence, we also neglect higher order 
contributions due to BIA. 


2. System with a Light-Hole Like Ground State 

According to Eqs. (12) and (13), on condition that 
6 > 2Ii^7r^72/(TOoL^) the ground state of the valence 
band is the first LH like subband. As in the case of a HH 
like ground state, we do not obtain a z-component in the 
effective SO field. However, in first order Lowdin pertur¬ 
bation theory, Eq. (B5), we obtain an additional linear 
term rj- and a cubic term £lt proportional to 643®". Fur¬ 
thermore, terms appear in third order which couple the 
electric field £z with the Dresselhaus term proportional 
to 5®4®“, App. B3. Since the SOC is a small correction, 
these terms are much smaller than the one not mixing 
both factors. Thus, according to the previous case of a 
HH like ground state, we have 

A’kin.LH = Ekin-, (29) 

14ff,LH = Keff,-(Aiql2) -I- (30) 

The effective vector field Jl due to Rashba and Dressel¬ 
haus SOC yields 


Olh = (31) 

with the additional term 

nb = blf^{k^kl,-klky,Oy. (32) 
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C. Summarized Results 

In summary, by developing an effective 2x2 model for 
a 2DHG we worked out the dominant contributions due 
to strain (/?, 0), Rashba (A_r) and Dresselhaus SOC ( 77 , 
Xd) to be considered in Eq. (16). The interplay between 
strain and Rashba or Dresselhaus SOC yields additional 
terms that are linear respectively linear as well as cu¬ 
bic in momentum. Thereby, we find that in contrast 
to the Rashba contribution the according Dresselhaus 
term in the confined system appears already in second 
order Lowdin perturbation theory. In respect of finding 
a conserved spin quantity we extracted the effective vec¬ 
tor fields J^hh for a HH like ground state and Olh for 
a LH like ground state in Eqs. (21) and (31). The fields 
cover a wide parameter space. In the next section, this 
will allow for identifying conserved spin quantities which 
do not require parameter-configurations which are diffi¬ 
cult to realize in real materials (e.g., 73 = 0 in Ref. 21, 
72 = —73 in Ref. 20). Thus, it facilitates the detection 
of long-lived spin states in experiments. 

III. CONSERVED SPIN QUANTITY 

Following the analysis of Ref. 7, our goal is to identify 
a conserved quantity E which is directly connected with 
/c-independent eigenspinors. The general ansatz is 

E = sol 2 x 2 + s • cr- (33) 

For this quantity to be conserved, it has to fulfill the 
relation [E,'Heff] = 0 which is true for 

n^Sz = ^ySz =0 A flySx — ^xSy = 0. (34) 

We are going to prove that one can find two solutions of 
this problem given by 


For the in-plane strain, the direction condition basi¬ 
cally requires symmetric normal strain components Cxx = 
eyy and a non-vanishing shear strain component exy ■ This 
situation can be generated by (110) uniaxial strain.We 
demonstrate this explicitely in App. D 2 for an experi¬ 
mental setup by use of a piezo crystal as done by Habib 
et at in Ref. 36. 

As we will see in the following, the constraint on the 
wave vector k of persistent spin states is crucial since 
it reveals that we found no conserved spin quantity for 
the whole fc-space but only for the averaged Fermi con¬ 
tour. However, this constraint is not surprising when we 
recall the case of persistent spin states in 2DEG. If the 
SO terms are linear in the wave vector, the condition 
for the existence of persistent spin states is fulfilled if 
the Rashba SOC coefficient is equal to the one for the 
linear Dresselhaus term. In this special case, the SO 
field is collinear in the whole A:-space. If the cubic Dres¬ 
selhaus term is included, we cannot find a quantity E 
which commutes with the Hamiltonian "H at every wave 
vector, though. Nevertheless, if the relations resulting 
from [E,H] = 0 are Fourier decomposed, similar to the 
procedure in Ref. 37, and only the lowest harmonics in 
the azimuthal angle is considered, one finds a condition 
for long-lived spin states. In contrast to the case without 
the cubic contribution, the found symmetry is, however, 
bound to an appropriate energy.^® This can also be un¬ 
derstood by studying the spin relaxation rates in diffu¬ 
sive n-type wires with Rashba and Dresselhaus SOC. One 
does not only find an additional spin-relaxation term due 
to the cubic Dresselhaus but also, the linear Dresselhaus 
coefficient is shifted.^® This shift depends on the Fermi 
energy. 

Next, as in the previous section, we consider separately 
the case of a HH like ground state and the LH like ground 
state. For the sake of simplicity, we apply the following 
replacements: 


Ef — 'y ( y ( A C^y)a/3Ck;3, (35) 

k,fc=fc||p a/3 


with ^ = ±, if either the strain is absent or its direction 
fulfills 


a I ^ 

0 = ± — = Y —. 

4 ^4 


(36) 


Here, creates a HH(LH) in the spin state a = ± for 
rrij = ±3/2 {rrij = ±1/2). We assume that the Fermi 
wave vector k\\p does not deviate much from rotational 
symmetry and thus can be replaced by its angular aver¬ 
age value fc||F = This situation holds for ma¬ 

terials close to axial symmetry, i.e., 72 = 73 , and a small 
strain amplitude /3. Therefore, we transform {kx,ky}^ 
into polar coordinates, fc||i?{cos(:/?), sin((/j)}^. Thus, if a 
hole, with a spin state given by { 1 , ± exp(^i 7 r/ 4)}^/-\/2 
and k = k\\p is injected into the two-dimensional system 
(including spin-independent scattering processes) its spin 
is not randomized. 



(37) 

II 
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to 

(38) 

/3 = Bkiip, 

(39) 


(40) 


and for the HH like state additionally 

V+=Vol2{kl)XB,+ - (41) 

In contrast to the discussed effect of the cubic Dressel¬ 
haus in an 2DEG, we will find persistent and not only 
long-lived spin states. 


A. Conserved Spin Qnantity in Case of a HH Like 
Gronnd State 

Making use of the definitions above and setting n = 
tthh for simplicity, we obtain for Ej the following equation 
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according to Eq. (34): 

0 = (42) 

= (cos((^)-^Sin((^)) 

• (?^ + C) + X [? + n(77o/K^ - 1)] 

+ [B^K'^nr + 2x (1 + r + < (r -1))] 

• cos(v?) sin((p)} . (43) 

This equation is fulfilled independently of the polar an¬ 
gle if if the ratio between Dresselhaus and Rashba SOC 
strength n and the strain strength factor B satisfy the 
relations 


42 


nif = C- 


2 (l-fir) 


2(i-r)-(ex)r«:2 


(44) 


More peculiar solutions occur for T 7 ^ 1 or in the pres¬ 
ence of strain. For certain parameter configurations, the 
SO field becomes collinear on the averaged Fermi con¬ 
tour. As a consequence, this SU(2) spin rotation sym¬ 
metry gives rise to a persistent spin helix.® For a start, 
by setting r = —l,x = ^ = lwe can recover the solu¬ 
tions presented by Sacksteder et al. in Ref. 20. In this 
case, we obtain n = 0, which is obvious since Dresselhaus 
SOC was not considered in Ref. 20, and 5 = 1, which 
is consistent with their results. It is remarkable that the 
presence of the linear Dresselhaus term, which was not 
considered either by Sacksteder et a/., does not alter the 
result. Note that here, too, the solution for the conserved 
spin quantity is bound to an averaged Fermi contour by 
Eq. (39). Moreover, recalling Eq. (34), one finds for these 
particular parameters a conserved spin quantity for every 
direction 9, given by 


and 


S = sin( 0 )(T 2 , — cos{6)ay, (51) 


„(±) _ . /^x(4-«;2)±W 

- V 2iv 

where 


(45) 


W = ^/K^ -f 8 (1 -fi r) 770 - 8 (1 -fi 2 r) + 16, ( 46 ) 

and are also solutions. 

If the Ck term in Eq. 23 can be neglected (full expres¬ 
sion can be found in App. C) real solutions for B^^^ are 
only found for n G where 


A+ = 
A_ = 


0,2-\/i + 2r-2^/(1+ r) (r-^M , ( 47 ) 

2^1+ 2r +2^(1+ r) (r-^M.oo) ( 48 ) 


if 3(Asg„(j.^)) = 0 which usually holds as typically F > 0 
and 770 < 0. 

In absence of strain, that is, (3 = 0, the fomulas above 
yield the following requirements on the ratio n and for 
the Fermi wave vector to get 


s 

0 

III 

s 

II 

1 + 

(49) 

; ; j 9+ 3 + T 

= V 2^,Ao,., r ^ 

(50) 


In this scenario, the SO field even vanishes at a specific 
value fco of the Fermi wave vector in the axially sym¬ 
metric case, i.e., F = 1. For this purpose, both linear 
and cubic Dresselhaus contributions are crucial. Hence, 
this solution was not existent in our previous publication 
Ref. 21. The axially symmetric case demands a vanish¬ 
ing Rashba contribution, Ah^+/Ad,-i- ~ 0. However, for 
the most semiconductors F ranges from 1 to 1.5. Thus, 
the cubic Dresselhaus SOC strength has to outweigh the 
Rashba SOC strength, i.e., |Ad,-i-| > |Ar,+ |. 


which generalizes the result from the previous publica¬ 
tion. Yet, we emphasize that the condition F = — 1 is 
rather unusual for most materials. Regarding the real¬ 
ization of persistent spin states in experiments, it is of 
interest to analyze whether the constrains allow for real¬ 
istic (i.e., typical for HI-V semiconductors) parameters. 
Thus, we present a concrete example for the conserved 
quantities Ej. At this, we assume 73 — 72 > 0, Eq. (Al), 
and 7 i > 0 to hold and choose, as an example, F = 1.2 
for the plots of and B^^^, Fig. 1. In order to draw 
a general picture and for simplicity we neglected the lin¬ 
ear Dresselhaus contribution, i.e., rjo —>■ 0. In the range 
A_|_, i.e., ^ = X thsA, comprises realistic values for n (e.g., 
K = 0.3 for a confinement width L = 100 A and a small 
Fermi vector = 0.01 A ) solutions are displayed in 

Fig. l(a,b). The solutions with a small B value, si 1 
in Fig. 1(b), are preferable since in this case the defor¬ 
mation of the Fermi contour is small. In addition, it is 
reasonable to assume that n > 1, since the Dresselhaus 
contribution is usually larger than the Rashba contribu¬ 
tion. The domain A_ (^ = —x) shown in Fig. l(c,d) 
is less realisic. Having large values of k implies a large 
width of the quantum well and a high Fermi energy which 
leads to populations in higher subbands where the model 
loses its validity. Also, for a large Fermi energy the Fermi 
contour is strongly deformed as can be understood by 
examining the term proportional to / 3 ^ 7 | of the kinetic 
energy, Eq. (19). A spherical approximation becomes in¬ 
appropriate. Moreover, if a strong strain is applied to 
the sample the appropriate model Hamiltonian needs to 
include also the coupling to the split-off band. 

For future devices like the spin field-effect-transistor it 
is not only of interest to find persistent spin states. In 
fact, samples are favorable where the injected particles 
undergo only a well defined spin-rotation. Thereby the 
initial spin state, with k being a good quantum number, 
is not necessarily an eigenstate. Here, well defined means 
that the rotation only depends on the distance between 
the injection and detection position. For n-type systems 

















FIG. 1. (Color online) Parameter configurations for and 
which yield the conserved spin quantity Ej for 6 = X^/4 
in case of a HH like ((a)-(d)) and LH like ((e)-(f)) ground 
state. The ratio between the Luttinger parameter 73 and 72 
is set to r = 1.2 and the linear Dresselhaus contribution is 
neglected, i.e., rjo —>■ 0. The domain for k is A+ for ^ = x s-nd 
A- else. The dashed lines indicate the according asymptotes 
at large width of the quantum well, (a) and (b) 

and for 5 = x = — 1 - IfC = X = l holds only the 

sign of n is inverted, (c) and Bi'^J, (d) ni~} and Bi~J 
for 5 = —X = 1. Interchanging of ^ and x only changes the 
sign of n. (e) and B^'^^ which yield the conserved spin 
quantity E_ for C = X = “1 ^nd (f) ^ = —x = —1. The case 
5 = 1 reverses only the sign of n. 


this condition was already analyzed in Ref. 7. Concerning 
a 2DHG as described in this paper, a spin-conserving 
condition which is valid for spin states with arbitrary 
wave vector cannot be found: The condition is limited to 
the averaged Fermi contour. For these states we can find 
an additional condition so that their precession depends 
only on the distance. At this, a necessary condition is 
an elastic scattering from impurities. The corresponding 
effective vector field in the case where the Eqs. (44,45) 
hold has the structure given by 


= {kx + ^ky)ip^^'>{ka;,ky) 




(52) 


In the case where (fca, + ^ky)ip^{kx: ky) depends linearly 
on ki {k\\p is a constant), the mentioned spin-rotation is 
only distant dependent. Here, one finds 

Lp^^\kx, ky) 


fcyi.72Afl,+ 

K2(±5;^^n;-H4r-K2) 


2(4±5xW)(i + r) 


- (2 ± 5xW -b 6r)«:2 + + 58(r2 _ i)^^^ . 

k\\p J 

(53) 


Thus, the special case where a well defined spin rotation 
occurs can only be found if F = 73/72 = ± 1 . 


B. Conserved Spin Quantity in Case of a LH Like 
Ground State 


Analogously, it is possible to find conserved spin quan¬ 
tities if the ground state is LH like. Yet, the structure of 
the SO field is more complex since it contains additional 
first order term due to Dresselhaus SOC, Eq. (32). As 
stated above, an [ 110 ] uniaxial compressive strain leads 
to (5 < 0 and therefore cannot be used. It is commonly 
known that for in-plane biaxial tensile stress LH ground 
state can be created, but in that case txy vanishes^^. 
Consequently, combined strain effects are necessary to 
generate the required condition. Nonetheless, we stress 
that we do not demand a strong in-plane strain ampli¬ 
tude for identifying a conserved spin quantity. In fact, 
an appropriate tensor component e^z is necessary. This 
component is encapsulated in the splitting 

Ahi,in = {kl) [{n^ - 1) 7i -b 2 (n^ -b l) 72] - < 5 , 

ZTTT-o 

(54) 


and thus in the SOC strength. For simplicity, we set n = 
riLH- Hence, the calculation of the conserved quantity is 
the same as before and valid as long as the deformation of 
the Fermi contour is not excessively strong. In this case, 
we find for the parameters and the relations 


,(±) = ^ 6(F-1) _ 

3xr (sg) ' k2 - 25 (3 (F + 1) + 2Q) 

and 


(55) 


b(±)= (48-3^2 + 12) ± 7 . 

Gx V ' 6 rK 2 ■■■ 

where we defined 


(56) 


V = 
Q = 


^Jm{Q + 3 r)^ + 24k:2 (Q + er - 3 ) -b 9«:4, 

Aii^hi 

|(^lbhlL^o)| 


(57) 

(58) 


The parameters and are plotted in Fig. l(e,f) 
for F = 1.2 and Q = 1 which is equivalent to an energy 
shift 6 = 2 |(Aii_hi|^_Q)|. We only find real solutions 

for and In contrast to the HH like ground 
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state for a realistic system with F > 1 we do not find a 
conserved spin quantity if strain is absent. 

In the last part of this section we apply the insights on 
the conserved spin quantity to a prominent semiconduc¬ 
tor. 


C. Example: p-doped InSb 


We choose p-doped InSb as an example to contrast the 
strained case yielding a conserved spin quantity with the 
one of a strainless sample. We assume a confinement in 
[ 001 ] direction with a depth of L = 100 A. To guarantee 
a low filling we set k\\p = 0.01 A“^. The used parameters 
here are listed in App. E. Further, we assume the addi¬ 
tional splitting due to strain between HH and LH like 
subbands to vanish, i5 = 0. Choosing an external electric 
field of 8z = 1.6 kV/cm and a [110] tensile strain direction 
{f-xy > 0), i.e., C = X = Ij allows for a persistent spin po¬ 
larization in [110] direction. Since we can assume a HH 
like ground state, we apply Eq. (44) and (45) and ob¬ 
tain the parameter for the in-plane strain strength to be 
b\ = 0.74 and the corresponding ratio between Dres- 

selhaus and Rashba SOC strength to be = —20.7. 

In Fig. 2(d) the resulting effective SO field is plotted 
and compared to the case where the [ 110 ] stress is ab¬ 
sent, Fig. 2(c). A stream plot. Fig. 2(a,b) shows that 
without strain. Fig. 2(a), the vector field vanishes ap¬ 
proximately at fco = 6.9 X 10“^ A“^ which is illustrated 
by the blue dotted circle. We see that even though the 
condition Eq. (49) on ng for the spin-preserving symme¬ 
try in the strainless case is not perfectly fulfilled, i.e., 
I > riQ = —29.6, the location where the field disap¬ 
pears is still well described by fcg in Eq. (50). Addition¬ 
ally, there is one source in the vector field at fey = 0 . 
Including strain, Fig. 2(b), gives rise to two additional 
sources that are centered at the crossing of the Fermi 
contour and the [110] axis. This can be understood by 
considering the factor {kx + ^ky) in Eq. (52). At these 
two sources vector field components are suppressed which 
are not collinear with the [110] direction. The Fermi con¬ 
tours split due to the SOC. Without strain, they are only 
slightly deformed as consequence of the band warping. If 
strain is present, the deviation of rotational symmetry 
of the contours is enhanced. The deformation is most 
intense in the [IIO] and [110] direction. To guide the 
viewer’s eye, we give in Fig. 2(e,f) a detailed picture of 
the SO field acting on the outer Fermi contour. In the 
case of strain the vectors lie parallel to the [ 110 ] direc¬ 
tion. As discussed above, the plots also show that the 
preserved spin quantity is limited to the averaged Fermi 
contour. The vector field regions which are noncollinear 
are strongly suppressed, though. This leads to a reduc¬ 
tion of spin relaxation even in the case of a general spin 
state injected into the 2DHG. 



FIG. 2. (Color online) Comparison of the SO field without 
applied strain (left column) and the case where strain gives 
rise to a conserved spin quantity E in case of a HH like ground 
state. In the presented case the SOC strengths are Ad,-i- = 
-42.5 eVA^ Ar,+ = 2.05 eVA® and ??+ = -34.1 x 10"®eVA. 
(a) and (b) Stream plot of the effective SO vector field. The 
green circle indicates the axially symmetric Fermi contour. 
The blue dotted circle corresponds to fco = 6.9 x 10“® A“^ 
where the field vanishes approximately, (c) and (d) Vector 
field and Fermi contours. The gray arrows indicate the spin 
polarization, (e) and (f) Detailed picture of the SO held that 
operates at the outer Fermi contour. 


Influence of linear Dresselhaus terms 


Moreover, we want to emphasize that in the cho¬ 
sen parameter regime the effect of the linear Dressel¬ 
haus contribution is only small. Fig. 3 shows how the 
field modifies if the linear Dresselhaus term proportional 
to r]+ is neglected. The strain-induced additional field 
sources move to a marginally lower Fermi wave vector 
by e = 3.1 x 10“^ A“^ where the conserved spin quan¬ 
tity is reobtained. To first order in 77 + the shift e can be 
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k, (A-‘) 

FIG. 3. (Color online) Stream plot of the effective SO vector 
field without linear Dresselhaus contribution, i.e., —>■ 0, in 
case of a HH like ground state. The strain-induced sources 
of the field move to a slightly lower Fermi wave vector by 
e = 3.1 X 10“® where the spin-preserving symmetry is 
recreated. 


£xx = £yy are essential. We have recovered the conserved 
spin quantity presented in Ref. 20 for the special case 
where 72/73 = — 1. In all circumstances, owing to the 
presence of both linear and cubic terms due to SOC the 
persistent spin states are bound to a Fermi contour. We 
have also demonstrated that only for this case and for 
72/73 = 1 one finds a spin rotation of a spin on the 
averaged Fermi contour which only dependents on the 
distance between the injection and detection position. 
Moreover, we have shown that for the existence of a con¬ 
served spin quantity in semiconductors which are acces¬ 
sible for experiments (e.g., systems with 72/73 « 1 ) the 
interplay between Dresselhaus SOC, Rashba SOC and 
possibly strain is crucial. 

In this way, shear strain has turned out to be a key 
component for an efficient manipulation of spin lifetime 
in 2D hole systems of zincblende structure. 
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Appendix A: Luttinger Parameter Relation 


We stress that the influence of the linear BIA terms be¬ 
comes even smaller for increasing Fermi wave vector kp 
and in other materials such as GaAs is less signihcant. 


IV. SUMMARY 

Summarizing, we identified conserved spin quantities 
in a (OOl)-confined two-dimensional hole gas in semicon¬ 
ductors with zincblende structure. Thereby, we derived 
the dominant contribution to the SO field due to Rashba 
SOC directly from an electric field which was miss¬ 
ing in our previous publication. Ref. 21. The significant 
effect due to Rashba SOC is only controlled by the sub¬ 
band gaps and not, as in the case of Dresselhaus SOC, 
by the conduction band gap. In view of recent publi¬ 
cations, we also included the effect of linear Dresselhaus 
SOC terms whose significance was pointed out in Refs. 27 
and 28 to be underestimated. The proper determination 
of the SO field enabled us to conclude that there are two 
possiblities for long-lived spin states. In respect of an un¬ 
strained sample such states exist only for heavy holes. It 
requires a certain ratio of cubic Rashba and Dresselhaus 
SOC strength defined solely by the Luttinger parame¬ 
ters 72 and 73 . Other spin-preserving symmetries oc¬ 
cur in presence of strain for both a HH like and LH like 
ground state. Here, a non-vanishing [110] shear strain 
component Cxy and a symmetric in-plane normal strain 


Here, we shortly focus on the Luttinger parameters 
72 and 73 . The warping of the valence band is directly 
proportional to their difference. Comparing with exper¬ 
imental results, for most semiconductors one finds the 
parameter 73 to be larger than 72 .^^^^"^ This becomes 
clear when describing the system using k • p method for 
band structure calculations, which yields the relation^^ 


73 - 72 


2 

3tooE(, 




(Al) 


Hereby, we follow the notation and phase conventions of 
Ref. 23 where Q S K and iQ being a momentum matrix 
element between the Fg^^^ valence and the r 7 c and Fgc 
conduction band states. The energy separation between 
the conduction band Ttc and the j = 3/2 valence bands 
is denoted as Eq > 0 .^^ 


Appendix B: Applied Approximations 


1. Lowdin’s Partitioning 

In this paper we start with a Hamiltonian which de¬ 
scribes HH and LH in the bulk. At the F point both types 
are degenerate. Imposing a confinement on the system 
reduces it to a quasi 2D system, generating HH like and 
LH like subbands. The simplification, which allows for 
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further analytical studies, is now to focus on the subspace 
spanned by either the set {|j = 3/2, mj = ±1/2) \n = 1)} 
or {|j = il2,mj = ±3/2) \n = 1)}, only, where n is the 
subband index. Let us call this subset A. Due to the 
confinement in growth direction and strain, this subset is 
well separated in energy at the F point from all other sub¬ 
bands (except the particular case where strain is exactly 
reversing the energy splitting due to the confinement be¬ 
tween the lowest band and the subsequent one at the 
r point). Assuming A having HH character, let us call 
the subset spanned by {\j = if2,mj = ±3/2) |(n > 1)), 
\j = 3/2, TOj = ±1/2) |n)} subset B. To end up with an 
effective model which focuses on the lowest HH(LH) like 
subbands, one treats the effect of subset B on the sub¬ 
set A in a perturbative manner. Here, one has to dis¬ 
tinguish between degenerate and non-degenerate pertur¬ 
bation theory since A and B may contain exact or ap¬ 
proximate degeneracies. Applying the quasi-degenerate 
perturbation theory, also called Lowdin’s partitioning, 
avoids this tricky task. It is described in great detail in 
the book by Bir and Pikus^° or by Winkler^^. With the 
subspace decomposition of our problem in mind, we give 
in the following a short formal description of the per¬ 
turbation procedure. Assume that T-L can be expressed 
as a sum of a Hamiltonian T-L^ with known eigenvalues 
En and eigenfunctions ji/n) and H' which is treated as 
a perturbation. Further, assume %' being a sum of a 
block diagonal matrix T-L^ with subsets A and B and Ti,^ 
describing the coupling of both subsystems, 

+ n' = n° + 71^ +n‘^. (bi) 

In accordance with Ref. 24 and 30, we define the indices 
corresponding to the states in set A, the in¬ 
dices l,l\l" to states in set B and the matrix elements 
between them as 


= (B2) 

Now one can find a non-block-diagonal anti-Hermitian 
matrix S which allows for a transformation of Ti yielding 
a block diagonal Hamiltonian % = . S and thus 

% can be found from a successive approximation. The 
first terms 


n = + -^( 2 ) -^(3) ^ 


are given by 


■^( 0 ) _ qjO 

^mm' ^mm' 5 

( 1 ) = 1^' 
^mm' ^mm' ’ 


nj (2) 
^mm' 


1 

2 




1 

Em — El 


1 

Em' — El 


-^(3) 

'^mm' 


2 ^ 

l.m" 


'IJ' 'IJf 'IJf 


{Em' — El) {Em" — El) 


± 


nL 


nj! nj! 


{Em - El) {Em" - El) _ 


(B3) 

(B4) 

(B5) 


(B 6 ) 



^[110] (A 


4x4 

14x14 


FIG. 4. (Color online) Absolute value |Ai?| of the BIA spin 
splitting of LH and HH states in GaAs for k| j [110] in the bulk 
system. The results are obtained by means of diagonalization 
of the full 14 X 14 matrix of the extended Kane model (or¬ 
ange) and the effective 4x4 model used in this paper (blue), 
following Ref. 24. Here, contributions from remote bands to 
the extended Kane model are excluded. 


2 ^ 

1,1' 


u' U' 'U' 


.{Em — El) {Em — Ell) 


{Em' - El) {Em' - Ei>) 


(B7) 


Note that each of the subsets A and B may be degen¬ 
erate but it is crucial that the subsets are chosen to be 
separated in energy, i.e., Em Ei. 


2. Dominant Invariants for the Cubic BIA Spin 
Splitting in Bulk Semiconductors 


Our starting point is the extended Kane model exclud¬ 
ing contributions from remote bands, i.e., Ck = 0. The 
remaining invariants for the Tg^, band block in this model, 
which give rise to cubic BIA spin splitting, are given by 
Eq. (3). In order to determine the coefficients, Lowdin’s 
partitioning in the energy gaps at the T point in the ex¬ 
tended Kane model can be applied. In Ref. 24 [Tab. 6 .3.], 
the coefficients in third order perturbation theory have 
been listed for several compounds. It reveals that in the 
bulk compared to 64 /®”, the terms proportional to 642 ®”, 
&®f" and 5®f" can be neglected. As an example, we com¬ 
pare the extended Kane model with the effective 4x4 
model used in this paper by calculating the absolute value 
|AA|. The latter is the BIA spin splitting calculated for 
both LH and HH states in GaAs for k||[110] in the bulk 
system. The result is plotted in Fig. 4 and shows very 
good agreement between both models. Deviations are 
only present at large k values. Note, however, that for 
a certain choice of parameters for the k • p model higher 
order corrections to the coefficients can be significant, as 
recently shown in Ref. 28. For a comparison of "Hbia 
with the terms used in Ref. 11 or Ref. 28 it is useful to 
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recast Eq. (3) for Cfc = 0 in the form 
^BIA = "(J • k) + ^ + (blf^ - blfn 

Oc 

1 


Vaka ( k, 


btfVHfn -1 


(B 8 ) 

with Vx = {Jx, Jy - Jl}, Ja and Kx = kxik^ - k^) (and 
corresponding terms). The relation between the coeffi¬ 
cients used in Ref. 11 and the one in Eq. (3) are thus 
given by 


^ 8 . 8 . ^ ^ 

0 


2i 


^ 8 . 8 . ^ _pp,Q 

= ‘^^pp'Q 


1 1 

n3 


5 

Eo ' 

Uo 

E'o 

+ Ag 

1 1 

r 1 


1 

Eo ' 

{e' 

E'o 

+ Ag 

1 1 

( 

1 


Eo ' 

{e' 

+ Ag 

E'o. 

1 1 

f 1 



Eo ' 

{e' 

E'o 

+ Ag 


13 , 


(B9) 


, ) — 


(BIO) 


(Bll) 


with 


and 



P = — (5|p,|A), 
mo 

(B13) 


P'= ^ {S\px\X'), 
mo 

(B14) 


Q=^)X\Py\Z') 

mo 

(B15) 

Ao = 

Q'? fi 

-^(^l[(VK))xp],|Z), 

(B16) 

A' = 

Q-? fi 

-^2myVo)xv\y\Z'), 

(B17) 

A- = 

-^(^l[(VK))xp],|Z'), 

(B18) 


3. LH like Valence Band Ground State: Mixing of 
the Electric Field and Dresselhaus Term 


Assuming a LH like valence band ground state and ap¬ 
plying Lowdin perturbation to third order, terms appear 
in third order which couple the electric field £z with the 
Dresselhaus term proportional to 64 )'®’', 


= - 

niix,rE 


— — 


= 0 . 

mix ,2 


256eg^&lf" 


(26: 


,8u8t>x /x2 
41 ^vK'^x 


256e£:^6tf" 


{2hlf'^ky{kl 


ky) + eSzkx) , 

ky) G-Szkx'^ J 

(B20) 


It is negligible if compared to those proportional to Ah^± , 
Ad,± or 77_, Eqs. (24)-(27), since we assume the SOC to 
be a small correction. Therefore, we neglect these terms 
in the calculation of the conserved spin quantity. 


Appendix C: Domain 


If the Ck term in Eq. 23 cannot be neglected the real 
solutions for are only found for k G £^sgn{^-x) '"^t^ere 


= 


A_ = 


0,2 J1 + 2r + (1 + r) - 2 ^( 1 + r) A 

V 2 Ad 


(Cl) 


2W1 + 2r + ^ ^(1 + T) + 2V(l + r) A, 00 


(C2) 


1 


with A = — (6C'^(1 -f T) -b 8y3CfcAD(l + 2r) (C3) 
2Af, 

+ MSkljziblf^ + 6 «f") + 32rAD)) (C4) 

and Ad = 72k^ipAD,+ - (C5) 


Appendix D: Strain 


with Vq the Coulomb potential of the atomic core, 
(X, Y, Z) the topmost bonding p-like valence band states 
and the antibonding s-like (S) and p-like (X',V',Z') 
states in the lowest conduction band. 

Notice that there is a difference in the application of 
perturbation theory to get the prefactors: According to 
Eq. [19a] in Ref. 42, e.g., one finds for Ag « 0 the term, 
which corresponds to 642 ®’', to be 


i,8vSv 

^42 


PP'Q 


1 

Eo{Eo + Eq) 


(B19) 


In contrast, Eq.(BlO) vanishes if Ag can be neglected. 


1. Deformation Potentials 

The deformation potentials a, 6 and d can be defined in 
different ways and there are several of them in the litera¬ 
ture. We list some relations between different definitions 
in Tab. I. 


2. Uniaxial Strain via Piezo Crystals 

Experimentally an unaxial strain can be conveniently 
implemented by the application of a piezo crystal since 
the deformation is tunable. In this setup, as done by 
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Eq. (4), Ref. 20 

Ref. 24 

Ref. 30,13,12’",23® 

a 

Dd - 

a+lb 

b 

- D 

-b 

d 

Id' 

-dfVS 


'^The sign of a has to be inverted. 
^The sign of b and d has to be inverted. 


TABLE I. Relations between different conventions for the de¬ 
formation potentials for the Fg^ valence band. 


Habib et al. in Ref. 36, the sample is fixed at one side 
of the piezo crystal where we align the poling direction 
of the piezo with the [110] direction. Depending on the 
polarity of the applied voltage, the piezo crystal extends 
(shrinks) along its poling direction and simultaneously 
shrinks (extends) perpendicular to it. Assuming the de¬ 
formation being completely transmitted to the sample, 
we can relate the strain coefficients of the sample, where 
the principal axes correspond to the three (100) axes, to 
the strain coefficients of the piezo, which can be directly 
measured. We define the transfered strain parallel to the 
poling and perpendicular to it as ej| and e(j_. Since due 
to Hooke’s law an in-plane strain generates also a finite 
out-of-plane component, the strain tensor of the sample 
becomes 


e = 


1 

2 


e,, - e I 


Cn - ej 

0 



(Dl) 


where C 12 and Cn are stiffness tensor components de¬ 
pending on the sample’s material. It becomes clear that 
in this situation the in-plane normal strain is symmetric, 
i.e., Cxx = iyy, and the shear strain component Cxy ^ 0 
as eJi and e(|_ have opposite sign, which corresponds to 
the situation demanded in Sec. III. 


Appendix E: BIA Parameters 


In Tab. II we list the coefficients for the invariants ap¬ 
pearing in the BIA Hamiltonian, Eq. (3), for some com¬ 
mon compounds (all values in eVA^, except for Ck in 
eVA)24. 



GaAs 

AlAs 

InSb 

Ck 

-0.0034 

0.0020 

-0.0082 

t8v8v 

O41 

-81.93 

-33.51 

-934.8 

t8v8v 

O 42 

1.47 

0.526 

41.73 

t8v8v 

O 51 

0.49 

0.175 

13.91 

t8v8v 

O 52 

-0.98 

-0.35 

-27.82 


TABLE 11. Expansion coefficients for the invariants in the 
used model, Eq. (3), up to third order in k which give rise 
to BIA spin splitting. All values in eVA^, except for Ck in 

eVA.24 
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